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OBJECTIVE 


The broad objective of this urban hydrology study of the Ross 
Creek Watershed was to assist the Santa Clara Valley Water District 
(SCVWD) with its reevaluation of their urbden hydrology 
1 
methodology . A key feature of the study was to be a new attempt to 
relate fitted hydrograph parameters to physical characteristics of 


the basin, and from this to develop a procedure to calculates design 


flows from the physical characteristics of the urban watersheds. 


More specifically, the major components of this study were to 
include : 

(1) Preparation of information and data on the drainage areas, 
precipitation gaging stations, stream gaging stations, 
precipitation, and streamflows; | 

(2) Characterization of the watershed, including its division into 
Subareas, pervious and impervious portions; 

(3) Investigation of the variation of routing (storage/discharge) 
relations with land use and size of drainage area; 

(4) Use of the Corps of Engineers computer program HEC-1 to 
determine, for each subarea, the nydrologic parameters which 
characterize Clarke's synthetic unit hydrographs for the per- 
vious and impervious portions; 

(5) Determination of regression equations which relate the above- 
mentioned hydrologic parameters to physical characteristics of 
the subareas, and from. those development of a procedure to 
calculate design flows from physical characteristics; 

(6) Verification of the design flow methodology on different 


storms; and 


(7) Writing a final report. 


EXTENT OF COMPLETION 


About half way through this project it became clear that step 
4 above, which tried to match complete hydrographs to observed 
ones, was far more complicated than originally anticipated. Near 
the three-quaters point, it became very apperent that the accuracy 
of the hydrologic parameters determined in step 4 was still not 
good enough to justify the step 5 regression study, and therefore 
thet the full scope could not be completed within budget. In a 
joint decision with the SCVWD staff, it was then decided to devote 
the remainder of the budget to improving the results of step 4 and 


to writing the final report (step 7). 


Some further improvement of step 4 was achieved, and steps 5 
and 6 were not attempted. Because full completion of step 3, and 
to a lesser extent step 2, depended somewhat on a good flow 
modeling capability which is not achieved in step 4, these two 
Steps were not completed to our satisfaction. For step 2, there 
only remains some uncertainty of the pervious and impervious 
allocations. .For step 3, it was only possible to investigate the 
sensitivity of the modeling process to changes in the 


storage/discharge relation. Steps 1 and 7 have been completed. 
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dy used the HEC-1 computer program of the Army Coros 
2 
of engineers to model urban hydrology behavior . Because cf tne 
complicated characteristics of urban areas, there are many parame 
ters and special situations which need to be considered in the 
calibration process. In our efforts to achieve the best simula- 
tions, we adjusted the parameters based on previous experience and 
on comparisons of hydrographs computed fron rainfall records at 
nearby stations with the observed hydrographs from flowrate records 


for the two gaging stations on Koss Creek. 


In the following sections, we briefly introduce the basin 
studied, its data, and the HEC-1 program. We first ran optimiza- 
tions on both subarea A (based on the flow records for station #21) 
and on the whole basin (based on flow records for station #51), to 
obtain guidlines for later adjustments of the key parameters 
including time of concentration, routing coefficient and loss rate 
parameters. . For the whole basin thos® parameters were optimized 
rather roughly because of the large area and hence more varied 
characteristics. For a smaller area such as subarea A, the optimi- 
zation can give better results, as discussed later. Following the 
optimizations, investigation procedures included (1) adjusting 
subarea boundaries, (2) testing alternative routing methods, (3) 
investigating the effects of various areal rainfall aieePi nations 
on the outflow hydrographs, and (4) testing and adjusting important 


controlling parameters, such as the three loss rate parameters, the 


cw 


hydrograph recession rats, the time of concentration, and the 


routing coefficient of Clarke's synthetic unit hydrograph, in order 
3 


to reconstitute the outflow hydrographs. 


The four procedures just mentioned 
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re further described and 
discussed in later sections. In appendicies this report also pre- 
Sants tables of parameters and constants used for the different 
approaches, and comparisons of computed and observed hydrographs. 


Information on eighty four computer runs is included. 


THE BASIN 


The Ross Creek basin is shown in Fig. 1. Its total area is 
9.96 square miles. The basin covers both a southwestern part of 
San Jose and the upper part of Los Gatos. It is divided into eight 
Subareas. There are rural, hilly areas in most of subarea A and in 
Small parts of subareas B and H. The other subareas, 1, J, L, M, 
and N, are all urbanized areas. Ross Creek originates in subarea A 
and discharges into the Guadalupe River. The basin studied, | 
however, ends at streamflow gaging station #51 in subarea N, about 
3,000 ft upstream from the Guadalupe River. Most of the contribu- 
ting urbanized areas are medium-density residential, except for 
small parts of subareas N ard M which. are high-censity residential 


areas. Some schools and parks are also included. 


The triangles on Fig. 1 locate the two stream-flow gaging 
stations, #21 and #51. The circles numbered 1 to 6 represent the 
flow collection points of the subareas; for example, point 3 is the 


collection point for subarea H and point 4 is the collection point 


for subareas I and J. 
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Fig. 1 The Ross Creek Walerbasin Map 


DATA & METHODOLOGY 


Input data are organized into three categories: precipi- 
tation data, parameters and constants, and streamflow data. 
Listing of the precipitation and streamflow data at 15 minute 


intervals were provided by the SCVWD. 


Precipitation data have been collected at several SCVWD rain- 
fall recording stations. From them, Stations 1, 36, and 125, whicn 
form 2 triangle covering the basin (see Fig. 1), were chosen to 
best provide local rainfall records at fifteen minute intervals. 
Records for ten complete storms were available. We chose more 
concentrated storms, which had clearer peaks, for our initial 
hydrograph parameter calibration purpose. Storm 1 was used 
initially in the calibration process to test, adjust, and then 
generate a good set of model parameters. Storms 2 and 5, which 
have similiar shapes to Storm 1, were chosen to check our adjusted 
parameters. The results are given in later sections. Fig. 2 shows 
the rainfall distributions for the ten storms based on basin ave- 


rage records. These distributions were used initially, with basin 
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pplied equally over the entire 
basin. Theissen and isohyetal rainfall distributions were later 
tested, in which the gage records for the different suboreas were 
weighted in accordance with the respective methods. The hydre- 
graphs generated from these two methods, as shown on the 

computer outputs, were actually only slightly different from eacn 
other but do improve the model parameters. from those obtained by 


using the basin average rainfall distribution. This ma 
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Table 1 ~-continued 
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comparing Runs 38 with 39, 40 with 41, through Runs 74 with 75. 
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hod was used in the final runs. Table 1 compares 


ct 


heissen me 
the rainfall amounts based on basin average, Tneissen and isonyetal 
methods for Etorms 1, 2, 4, 5 and 7. The area percentages covered 

by each of the three rainfall gaging stations for both the Theissen 


and isohyetal metnods are given in Table 2. 


The many parameters required for HEC-1 ara listed in Appendix 
A. Starting with Run 9, th? seven subaress EB through NM were 


divided into pervious and impervious parts. This approach allowed 
different parameters values to be prescribed for the pervious and 
impervious parts of gach subarea, providing more accurate 


simulations. 


Streamflow data for both Stations #21 and #51 were read in 
only to enable the "observed" hydrographs to be plotted and com- 


pared with the computed ones. 


The hydrograph for the whole basin was generated by the 
following procedure: (1) generate the hydrographs for the pervious 
part and impervious part of the most upstream subdarea, (2) add the 
two hydrographs and then route (by the Modified Puls method, des- 
cribed later) the result through the subérea, to form tne cutflow 
hydrogragh for this specific subarea, (3) route the hydrograph 
along the creek channel using the Muskingum method (discussed in a 
later section) to reach the catch point, and (4) if there are ary 
ether hydrograghs also reaching that catch point, add them together 


to form the multi-subarea hydrograph for the total area upstream of 
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that catch point. RKepeat these four. procedures until the ia 
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Table 2 THE SUBAREA CHARACTERISTICS 
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caten point is reached. Since this is the position of Gaging 


Station #51 and the cutlet point of the wnele basin, we thus obtain 


the outflow hydrograph for the whole basin. 


There are eighty four selected computer outputs in Appendix B 
which display the hydrcegraph comparisons and their data for differ- 


ent approaches. Each "run" was numbered for easy later reference 


COMPUTER PROGRAM HEC~-1 


The HEC-1 model deveiopsd by the Hyrologie Engineering Center 
of the Army Corps of Engineers is designed to simulate the surface 
runoff response of a river basin to precipitation by representing 
the. basin as an interconnected system of hydrologic and hydraulic 
components. The HEC-1 program used was tne 1931 version, which has 
completely revised input and output formats (see the 1981 HEC-1 
User's Manual ). HEC-1 has some limitations, the mejor one being. 
that simulation is limited to a single storm due to the fact that 
provision is not made for soil moisture recovery during periods of 
no precipitation. The storms tested in this study, presumably 
because of this limitation. and the possibly alresady-saturated soil, 


yielded less than satisfactory predictions. 


No modification was made to the HEC-1 program other than to 
the input/output code, which was changed to suit Santa Clara 
University's DEC-2060 system. Three program functions were 


employed in this study: (1) unit hydrograph and loss rate parameter 


optimization, (2) watershed unit nvdregraph generation, and (23) 
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comparison between computed and observed hydrographs. 


appendices of this report. 


included in Appendix B. 


OPTIMIZATION 


The optimization function cf HEC-1 provides a set of opti- 
mized: parameters based on observed rainfall and runoff records. 
The parameters optimized in the study were the time of concentra- 
tion Te, routing coefficient R and the loss rate parameters STRKR, 
DLTKR, and RTIOL. The optimization results are shown in Runs 1 and 


2 of Appendix 8B. 


Because there are eight subareas in the whole basin, the 
optimization results yield "lumped" or averaged parameters which 
can only serve as guideline for calibrating the subares parameters. 
Although the error of the sum of flows and the average absolute 
percentage error were low, this doesn't mean that we can apply the 


optimized paérameters directly to represent all the subar 


om 
ie) 
177) 


Adjustment based on our knowledge of the characteristics of each 
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subarea was required in order to achieve more realistic model 
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parameters. This adjustment procedure is discussed late 


Parameter Adjustment and Reconstitution. 
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SUBAREA DELINEATION AND HYDROGRAPH GENERATION 


Tne reason for dividing the watersned into suberesas was to 


the basin's physical 
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characteristics. Eacn subarea should be an individual area with no 
flows across bounderies except at the outlet. For the Ross Creek 
watershed, larger storms sometimes fill up the stormdrain system 
and runoff then travels in the gutters, sometimes in different 
directicns from the stor drains. As a result, the subarea 


boundaries are changed for an uncertain time and for an uncertain 


Fraction of the runoff. These dynamic or moving boundaries, 


ct 


although occurring, were not treated as an important factor. This 
phenomenon may explain why simulations sometimes did not closely 


match the observed record. 


Stormdrain maps from the City of San Jose and the Town of Los 
Gatos were examined. We found that the original Subarea map pro- 
vided. by the SCVWD (see Fig. 1) had some bounderies which did not 
agree with the subareas defined by the storm drains. Early studies 


of ways to model drainage area characteristics included the 


183) 


following two approaches: (1) based on.stormdrain system maps, 
further divide the subareas so that no storm water flows across the 
boundaries, and (2) a micro-area approach, which treats each home 
lot as an individual area, and in which other land uses such as 


parks and schocls also generate individual hydrographs. Ideally, 


all the individual hydrographs in the basin can de added up to 
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We nad to abanden our efforts on these two tedious approaches 
because of their heavy time and manpower needs. However, Approach 
(1) is recommended for further study in the Ross Creek aree decnuse 
there have been some changes such as blocking off pipes and new 
construction in the area during the last decade, although we found 
the original map (Fig. 1) to be generally correct. The only change 


we made was to eliminate the srea downstream of gaging station #51. 


h. 


The subareas were divided into pervious and impervious perits, 
because? parameters such as Te, R, the less rate parameters, and 
consvancs such as baseflow aoa area differ considerably between 
these two components of a single subarea, but are similiar to their 
respective parts of other subareas. This provides a2 more 
reasonable approach to hydrograph generation. The representation 
of separate pervious/impervious parts also has the advantages that 
we cen describe the subarea more precisely and we can clearly 
define the impervious fraction of a specific subarea. Hote that 
‘the parameter RTIMP, ratio of imperviousness, should not be used 
under such @n approach however, or this will result in a double- 


accounting for impervious areas. 


ROUTING METHODS 


Flows in the main creek channel were routed using the Muskin- 
gum method. The creek was divided into five reaches, defined by 
collection points 1 through $ (where point 5 was taken to be at 


Station #51) shown on Fi 1. Thus the flows which entered tne 


oa 


upper coliestion point were routed through the reset co the lover 


collection point, where inflows from other sources could be added 
in to obtain the input for the next reach. The veius of 
Muskingum's K is approximately equel to the flow travel time 
through the reach. The average flow velocity used to allcecats K 
véelues was 4.5 fps. Muskingum's X is 2 weighting factor lyin 


between O and 1, and the following relationship between K and X 


should be satisfied: 


The X value was assumed to be 9.1 as recommended in the previous 
‘Study . A value of 0.3 was tested (see Run 21), but this chan 
was found to have negligible effect on the hydrograph. 

Flows through the stormdrain systems were routed using the 
Modified Puls routing method. It detained the water volumes and so 


tne rodified 


flattened the peaks of the hydrographs. We applied 

Puls routing to the whole basin, to subareas B through N, and to 

‘Sudareas I through N only (see Runs 3 to 8). The shapes of the 

resulting hydrographs showed that the larger the area covered by 

Modified Puls routing, the flatter was the hydrograph generated. 
The storage-discharge relationship used for the Modified Puls 
; 1 


routing method was taken from previous work . This relationship 


was changed several times to experiment on its impact on the hydro- 


graphs. Eighty percent and fifty percent of the storage in tne 


3 2 


the water volumes stored as expected, but the original reletionsnio 
? > By 


was used in the latest runs because it produced a more accurate 


hydrograph. However, this relationship should be «adjusted to 


reflect the updated storage-discharge curves presently being deve- 


leped by SCVWD. This we felt will directly affect the hvirograph 
which is worth exploring, because of the sensitive reaction from 


the different relationships. More study of tnis is recommended. 


RAINFALL DISTRIBUTION 


The distribution of earnratd within a storm was found to 
affect the hydrograph only slightly. Three types of precipitatio 
distribution have been tested. Basin average rainfall used the 
same rainfall over the whole basin, while the Theissen and 
Isonyetal methods yieided different rainfalls over each subere2a. 


For storms #1 and #2, the subarea Isonyetal method gave the best 


prediction, followed by the subarea Theissen method and lastly th 


’ 
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basin average rainfall. There was a difference of about 2% between 


ct 


the basin average and the Theissen methods, and about 12% 
difference between the Isohyetal and the Theissen methods in 
average absolute percentage errors. Both the Theissen and the 


Isohyetal methods were tested for Storms #4, #5 and #7 in Runs 38 


through 42. Table 3 shows the difference between Theissen and 


#4, #5, and #7. | . 
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isohyetal rainfall distributions on subareas for the storms #1, #2, 


Table 3 — accumuLatep RAINFALL OF STORMS #1, #2, #4, #5 and #7 
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Table 4 — THE OBSERVED LOSS RATE OF BASIN AVERAGE RAINFALL 
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PARAMETER ADJUSTMENT AND RECONSTITUTION 


The optimizstion procedur? decribed sarlier yielded "lumped" 
or average parameters for the entire basin. From our knowledge of 


elear that most 
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the physical charact2ristics of the basin, it + 
of these parameters should vary from area to area within the basin. 
The objective of parameter adjustment was to identify for each 
subarea parameter values, usually different, which were more reali- 
stic and which together deneruied improved, reconstituted outflow 
hydrographs. This process of identifying improved parameter values 


is referred to as calibration. 


Subarea parameters which had the most significant effects on 
reducing the error in the calibration were found. to be (1) RTIOR, 
which defined the recession rate of the subarea baseflow, (2) 
STRKR, which defined the starting loss rate of a subarea, and (3) 


DLTKR, which defined the accumulated loss while the loss rate kept 


increasing. 


RTIOR adjusted the shape of the recession curves which 
enabled the calculated curves to match well with the observed curve 
(Run 17). But for a different storm, it produced too fast a 
recession (see Run 40). STRKR and DLTKR control the total 
hydrograph volume and affect the accuracy of the hydrograph repro- 
duction a lot. According to the characteristics of the impervious 
and pervious surface, we have used higher loss rates parameters for 
the impervious parzv than the pervious part. This approach helped 
the computed hydrograph curves match the observed curves fairly 


well. In Run 22, we further reduced both parameters for the per- 
30 


vious part to very small amounts, and then the second nydrograph 
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peak appeared. This suggests that the sensitivity of the simula- 
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tion response improved. 


The other paremeters changed were Te and R of Clarke's uni 
hydrograph. A set cf Te and R values was obtained by trial and 
error. We then odtained a second set dy hyrograph analysis of the 
observed hydrographs for storms #1 end #2, and Saw mo significant 
improvements in hydrograph reproduction (Runs 34, 34-1, 35 and 35- 
1). We also conducted a sensitivity test for the routing coeffi- 
cient ’R with a 300% change and obtained undesirable results (Runs 
11 and 12). This suggests that we don't expect Te or R to have 
beneficial influences and that the values used for them are in 


reasonable ranges. 


The parameter adjustment tests were made with five storms: 
#1, #2, #4, #5 and #7. These represented different rainfall and 
runoff situatuions. Storms #9 and #10 were impossible LO test: 
Properly because they lacked rainfall records at Staticn 1, which 
affected the Theissen and isohvetal rainfall distributions over the 
subareas. Storms #1 and #2 are similiar, so the parameters were 


calibrated around these two storms, for which the finel accuracies 


of reproduction. were (see Run 32 and 33): 


Storm #1 Storm #2 
Sum of Flows: ; -3.07% ~6.17% 
Standard Error: - 22u 150 
Ave. % Abs. Error: U9.97% 42,454 


After developing the best set of parameters, we tested them un- 


changed on Storms #4, #5 and #7. Storm #4 had a fairly even dis- 
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tribution of rainfall, In Storm #7 the rainfall intensities rose to 
a sharp peck, and Storm #5 fell between the extremes of Storms #7 


and #4. The results (Runs 39, 41 and 43) were: 


Storm #4 Storm #5 Storn #7 
Sum of Flows: THY. Cus -0.79% 116.605 
Standard Error: 26? 125 210 
Ave. % Abs. Error: 236.673 150.625 262.65% 


These results Suggest that, at this point, we cen only satis- 
factorily apply the parameter set to storms similiar to Storm #1 
with short, intense periods of rainfall. Storm #5 showed a better 
reproduction than both Storms #4 and #7 becsuse it is more similiar 
to Storms #1 and #2 than the other two storms. Storms #3, #4 and 
#8, belonging to the "distributed" rainfall group, should be calibd- 
rated end tested together separately. It seems that we might need 
two or even more different parameter sets for the basin for the 
different storm types. In runs starting with Run 44, the 
adjustment and improvement were made mostly by changing hat ances 
STRKR, DLTKR and Te. It was evident that we should match the 
computed and observed total volumes defore we considered the shape 
differences. In Runs 44 to 53, we adjusted STRKR and DLTKR by 
increasing them slightly, and expected a lower total volume on the 
computed hydro-graph because of the higher loss rate. The results, 
as expected, gave lower average absolute errors. | DLTKR, 2s the 

al described, should be storm dependent. In later runs (Runs 
76 to 84), it was set at a higher value to obtain better 
comparisons between the equivalent depth of tne computed and 


observed hydrographs (see Table 4) because there was still 


insufficient loss indicated. We used this comparison as 2 
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guideline to adjust tne loss rate parameters. 


DLTKR from storm independent to storm dependent in later 


did not help much to reduce the errors. 


ie) 


Although we changed 


runs, this 


Dy 


ecmparing Runs 44 with 54, 


with 62. 


Reducing the delay in time to 


We obtained a low volume difference 


puted hydrographs. This delay time 


problem. We were never able to get 


respond fast 


hydrograpgns. 


the range of 0.75 to 1.5 hours, depending on the storm. 


and 37 we reduced STRKR for 


faster rising curves 
value of STRKR even 


reduced the 


loss. 


Sate for the overall 
tation that the reduced initial loss 


45 with 55, 


enough to match the rising 


each subares 
but without success. 
more 


This adjustment was based on 


46 with 56, through Runs 52 


The results were almost the same. 


peak was always the goal after 


between the observed and -comn- 


has always been a significant 


the computed hydrographs to 


curve of the observed 


The peak lag caused by this response delay was in 


In Runs 236 


in the hopes of achieving 


In Runs 64 69 we 


inerease 


Ou 
o 
rm 
= 
- 
oe] 


and TO compen- 


the expec- 


rate would 


cause the computed 


hydrograph to rise faster and that we could meawhile maintain the 


total volume ebout the 


Reducing STRKR by 50% and inereasing DLTKR. by 100% 


average absolute error to increase 
Theissen rainfall distributions of 


by about 50% for subarea isohyetal 


error on Storm #5. 
less err n Storm #5 


still delayed without 


In Runs 70 to 75, we reduced 
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same by adjus 


ting 
Dy about 


storms 


Fainfall, distributions, 


Te 


DLTKR simultaneously. 
caused the 
100% for the subarea 
#1 to increass 


and #2 and 


Dut gave 


The computed hydrograph peak time was 


any significant change. 


by 50% and expected to see 


some peax shifting. Although 


Te was then very smell, this shifted 


the peak by only about 0.25 hr. = 15 minutes. 


y 


In order to check the Te obtained oy trial and error, we then 
used the Kirby-Hathaway equation to calculate the Te for each 
Subarea by specifying the average travelling length and the 
Manning's n -values and slopes as shown in Table 5. In addition, 
as mentioned before, PLTKR also was increased in Runs 76 to 384. 

The results showed slight improvement. Runs 77, 80 and ¢3 compared 
computed and observed hydrographs for subarea A only. nese 
Suggested that our estimation of loss rates might be too high, 
because the volumes are significantly less than the observed 
volumes for subarea A. We then used the observed hydrograph at 
Gaging Station #21 anc wonvinued the modelling downstream to Sta- 
tion #51. The results (Runs 78, &1 end @4) showed good improve- 
ments in both the time to the peak and in the volume of the 


computed hydrograph (30% increase). 


- CONCLUSIONS 


Differences between the volumes of cbserved and computed 
hydrgraphs may actually be legitimate because of some doubtful 
flowrate records at Station #51. (These doubtful records were 
noticed at the end of this study when comparing flows at Stations 
#21 and #51.) For example, in Run 7@ we have a percentage error on 


the sum of flows as 37.21%, in Run 81 it is 32.61%, and Run 84 has 


Tc,= Kirby-Hathaway; 
Table 5. THE Tc COMPARISON, where: Tc.= Graphic Analysis; 
Tcs= previous best set from trial & error 


(ee ee eee oe 
Ge pee oe mee oe 
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we adjust other parameters. We also conclude that the stubborn 


~~ 


nature of the peak lag problem Suggests that the observed data may 


not have been synchronous. 


We hope that the efforts made on this study and the 


discussions held will be helpful for future research and 
development in the field of urban hydrology simulation. It is 


a 


clearly avery complicated system and there s 


ct 


ill are many problems 


needing further study. 


RECOMMENDATION 


The testing and adjusting procedures performed and results 
cbtained in this study lead to the following pe sounenaseion Ss 

(1) Manning's n needs more careful adjustment because it showed 
an important effect on Getermining Te; 

(2) parameter DLTKR plays an important part both in volume adjust- 
ment and in the average absolute errors. Because there is no 
clear definition for DLTKR and a lack of guidelines, more study 
is needed in order to find a better approach to adjust it 

(3) the storage-discharge relationship directly controls the volume 
and the shape of a2 Ayaromrach® Some improvements to the old 
data should be made to reflect subsequent changes in the 
stormdrain system in the urban. area, which should improve 


simulation. 


< 


(4) further divide the subarea so that there are no storm flows and 


drainage accross the subarea boundaries. 
38 
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